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Resumen 


CAMPOS-ARANDA, D.F Obtención de 
empíricas para estimación de crecientes de diseño en la 


ecuaciones 


Región Hidrológica Núm. 10 (Sinaloa) de México. Tecnología 
y Ciencias del Agua. Vol. V, núm. 1, enero-febrero de 2013, pp. 
125-143. 


Las crecientes de diseño (Q,,) son la estimación hidrológica 
más frecuente debido a que son básicas para dimensionar 
y revisar todo tipo de obras hidráulicas. Su estimación en 
cuencas sin datos de gasto máximo anual se basa en los 
métodos regionales, que emplean tal información de una 
región o zona geográfica para obtener resultados que son 
aplicables dentro de ella. Los procedimientos regionales 
son de dos tipos: los similares al índice de crecientes y el 
desarrollo de ecuaciones de empíricas. Los primeros se basan 
y deben su capacidad de predicción a la homogeneidad de la 
región; en cambio, el segundo método permite cierto grado 
de heterogeneidad, la cual es tomada en consideración en los 
propios modelos de regresión generados. Empleando como 
variables predictivas el área de cuenca (4), la longitud del 
cauce principal (Lcp) y las predicciones de lluvia máxima 
diaria (PMD), se establecieron cuatro ecuaciones empíricas, 
seleccionando dos, según sus índices de desempeño, para 
cada uno de los seis periodos de retorno establecidos. Se 
procesaron 23 registros de gastos máximos para obtener 
los valores de Q,, con base en la distribución general de 
valores extremos. Para la obtención de los coeficientes 
de ajuste de cada ecuación se emplearon tres métodos: 
el de mínimos cuadrados de los residuos en el dominio 
logarítmico, su corrección por sesgo y el de optimización 
numérica minimizando una función objetivo. Se obtuvieron 
seis ecuaciones empíricas para cada uno de los seis periodos 
de retorno analizados. El contraste numérico en seis 
estaciones hidrométricas no utilizadas en la deducción de las 
ecuaciones y en seis de las procesadas, las de mayor periodo 
de registro, indica una regularidad en las estimaciones y 
destaca la capacidad predictiva de las ecuaciones empíricas 
desarrolladas. Se recomienda aplicar el procedimiento 
propuesto en otras regiones para encontrar este tipo de 
ecuaciones, que permiten estimar las crecientes de diseño de 
manera sumamente fácil. 


Palabras clave: regresión lineal múltiple, error relativo, 
índices de desempeño, algoritmo de Rosenbrock. 


Abstract 


CAMPOS-ARANDA, D.F Obtainment of empirical equations for 
design floods estimation in the Hydrological Region No. 10 (Sinaloa) 
of Mexico. Water Technology and Sciences (in Spanish). Vol. V, 
No. 1, January-February, 2013, pp. 125-143. 


Floods designing (Q,,) is the most common hydrological estimation, 
because it is basic to sizing and reviewing all types of hydraulic 
works. Its estimation in watersheds without annual maximum 
flow data is based on regional methods, which use such information 
from a region or geographical area, in order to obtain results that 
apply within it. There are two types of regional procedures: those 
related to the index flood and those that develop empirical equations. 
The first type is based on the homogeneity of the region, hence its 
predictability; whereas the second method allows some degree of 
heterogeneity, which is taken into account in the regression models 
generated. Four empirical equations were established, using as 
prediction variables: the watershed area (A), the main channel length 
(Lcp) and predictions of maximum daily rainfall (PMD). Only two 
equations were selected for each of the six return periods established. 
23 records of maximum annual flow were processed to obtain 
values of Q,, based on the General Extreme Values distribution. To 
obtain the fit coefficients of each equation three methods were used: 
the least-squared residuals in the logarithmic domain, their bias 
correction and the numerical optimization minimizing an objective 
function. Six empirical equations for each of the six return periods 
analyzed were obtained. The regularity of the estimations was tested 
by numerically contrasting results from six gauging stations (not 
used in the derivation of the equations) and six of those processed 
(those with the longest period of record). Such analysis allowed to 
point out the predictive ability of the empirical equations developed. 
It is highly recommended to apply the proposed procedure in other 
regions to find such empirical equations for design flood estimation, 
because of their very simple application. 


Keywords: multiple linear regression, relative error, performance 
indices, Rosenbrock algorithm. 
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Introducción 


Las crecientes de diseño, conocidas de manera 
coloquial como avenidas máximas, son la esti- 
mación hidrológica más común, debido a que 
con base en ellas se dimensionan, en las etapas 
de planeación y/o diseño, todas las obras de 
infraestructura hidráulica; además, permiten 
la revisión de la seguridad asociada con tales 
Obras cuando ya están construidas, sean éstas 
de aprovechamiento o de control (embalses), de 
protección (diques y encauzamientos) o de cruce 
(alcantarillas y puentes). El procedimiento más 
confiable para la estimación de las crecientes 
de diseño consiste en el ajuste de un modelo 
probabilístico al registro disponible de gastos 
máximos anuales, para obtener estimaciones o 
predicciones relativas a diversas probabilidades 
de excedencia, cuyo recíproco es el periodo 
de retorno o intervalo promedio en años entre 
la ocurrencia de un evento igual o mayor. El 
procedimiento citado se denomina de manera 
genérica análisis de frecuencia de crecientes 
(AFC) y actualmente emplea como modelos 
probabilísticos prescritos la distribución Log- 
Pearson tipo III, la General de Valores Extremos 
y la Logística Generalizada, ajustando éstas a 
través del método de momentos, de máxima 
verosimilitud o de momentos L. 

La desventaja fundamental del AFC es la 
escasa disponibilidad de registros de crecientes 
en los sitios de interés. Por ello, se recurre al 
análisis regional de frecuencia de crecientes 
(AREC), el cual utiliza la información de gastos 
máximos anuales disponible en una región 
o zona que puede ser considerada homogénea 
desde un punto de vista hidrológico. Las 
técnicas de ARFC utilizan de manera conjunta 
la información citada y sus resultados resultan 
aplicables a cualquier sitio o localidad ubicada 
dentro de la región homogénea. Por ello, el 
ARFC permite estimar crecientes de diseño en 
cuencas sin aforos. Entre los métodos clásicos 
del ARFC se tienen la técnica de las estaciones- 
años, el índice de crecientes, el ajuste de 
modelos probabilísticos con base en parámetros 
estadísticos regionales y la obtención de 


ecuaciones de regresión potencial, que rela- 
cionan las crecientes de diseño obtenidas en 
cuencas aforadas, con ciertas características 
fisiográficas o climáticas de éstas (Cunnane, 
1988; GREHYS, 19%6a; Escalante-Sandoval y 
Reyes-Chávez, 2002). 

El GREHYS (1996a) indica que se puede 
establecer una distinción clara dentro de las 
técnicas del ARFC; por una parte, se tienen los 
métodos similares al índice de crecientes y por 
la otra al método que desarrolla ecuaciones 
de regresión empíricas. Los primeros se 
basan y deben su capacidad de predicción 
a la homogeneidad de la región y por ello la 
definición más precisa posible de ésta resulta 
fundamental. En cambio, el segundo método 
permite cierto grado de no homogeneidad, 
la cual es tomada en consideración en los 
propios modelos de regresión desarrollados. 
Entonces, si la región bajo estudio puede 
ser caracterizada con un alto grado de 
homogeneidad, los métodos semejantes al 
índice de crecientes serán preferibles; por el 
contrario, si existe incertidumbre al definir las 
regiones homogéneas o la zona estudiada está 
conformada por varias subregiones y no existe 
suficiente información hidrométrica en ellas 
para ser trabajadas de manera independiente, 
el método de las ecuaciones empíricas será la 
mejor opción dentro del ARFC. 

El objetivo de este trabajo consistió 
en desarrollar el método regional de las 
ecuaciones de regresión empíricas en la 
Región Hidrológica Núm. 10 (Sinaloa) de 
México, para la cual se ha encontrado que 
está integrada por dos subregiones: una al 
norte y otra al sur, según Escalante-Sandoval 
(1998); o bien por tres subregiones: zona norte, 
de montaña y planicie costera (Gutiérrez- 
López et al., 2004; Campos-Aranda, 2008a). 
Se procesaron 23 registros de gastos máximos 
anuales, con amplitudes que oscilan entre los 
19 y 56 años, con 36 años como valor medio. 
Se desarrollaron cuatro tipos de ecuaciones 
empíricas, cada una para seis periodos de 
retorno y todas se generaron con dos métodos 
básicos de ajuste, los mínimos cuadrados 
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de los residuos en el dominio logarítmico y 
la optimización numérica multivariada no 
restringida, aplicada a través del algoritmo de 
Rosenbrock. Como las ecuaciones obtenidas 
por mínimos cuadrados fueron corregidas 
por sesgo, se generó un total de 72 ecuaciones 
empíricas, de entre las cuales se seleccionaron 
36, por tener sus mejores indicadores de 
desempeño, éstas se probaron en seis 
estaciones hidrométricas no utilizadas en su 
deducción y en seis de las procesadas, las de 
registro amplio. Las conclusiones formuladas 
destacan la regularidad de las predicciones 
estimadas y sugieren el desarrollo de este 
tipo de modelos empíricos en otras regiones 
del país con base en los tres planteamientos 
de ajuste seguidos. Tales ecuaciones empíricas 
permitirán estimar fácilmente las crecientes 


de diseño en otras zonas geográficas del país. 
Desarrollo 
Ecuaciones de regresión potencial y su ajuste 


En la ingeniería hidrológica es común el uso 
de ecuaciones empíricas para realizar estima- 
ciones en sitios donde no existe información 
hidrométrica. En relación con la estimación de 
crecientes de diseño, la regresión potencial ha 
sido el modelo más utilizado para describir la 
relación entre las predicciones (Q,,), y ciertas 
características fisiográficas y climáticas (X) 
de las cuencas de las cuales proceden las 
estimaciones involucradas. Entonces, para una 
cierta región considerada homogénea, se tiene 
(McCuen et al., 1990; GREHYS, 1996a): 


Qr, =p X>-X2-.Xóm-8 (1) 


en donde By Bj 8 
ajuste y e es el error multiplicativo del modelo. 


son los parámetros de 


Para estimar los parámetros b, ha sido práctica 
común linealizar la ecuación (1), tomando 
logaritmos y entonces estimar los coeficientes 
de una ecuación de regresión lineal múltiple 
(RLM): 


log Q7, =l0gb, +b, :logX; + b, :logX, 
+-*+b, log X,, +log e (2) 


haciendo y = log Q,,, 1, =1l0g b,, a,=b, y x,=1l0g8 
X, con ¡= 1, 2,..., m, se obtiene la ecuación de 
la RLM: 


y =0)+b,:x,+b,:X +0 +b,,X,, + (3) 


cuya solución matricial de mínimos cuadrados 
de los residuos (MCR) se puede encontrar 
en cualquier texto de métodos numéricos, 
por ejemplo, en Campos-Aranda (2003). El 
coeficiente b, es igual a 10% o bien a e”, según si 
se utilizaron logaritmos decimales o naturales; 
los demás coeficientes tienen correspondencia 
directa. Cuando el error en la ecuación (1) es 
de adición, la RLM no se puede utilizar para 
estimar los parámetros de ajuste y entonces la 
optimización numérica ha sido empleada. 


Corrección por sesgo para el ajuste de MCR 


Desde principios de la década de los años 
setenta se comenzaron a desarrollar métodos 
para corregir los resultados del ajuste de la 
ecuación (1) a través de una ecuación de RLM 
debido a que el ajuste de mínimos cuadrados 
de los residuos minimiza a éstos en el dominio 
logarítmico, al emplear los logaritmos de las 
variables predictivas (X,) y de la dependiente 
(Q,,), como se estableció en la ecuación (2). 
Lo anterior implica que el coeficiente de 
determinación (R?) encontrado está evaluado 
en el dominio logarítmico y por ello se cita 
después como log R?. 

McCuen et al. (1990) exponen y aplican tres 
procedimientos para corregir el coeficiente 
b, también llamado ordenada al origen, por la 
similitud entre la ecuación (3) y la fórmula de la 
línea recta. El primer método corrige la mayor 
parte del sesgo y su planteamiento consiste 
en obtener una estimación de b, no sesgada, 
haciendo cero los residuos (e) en el dominio 
real, es decir: 
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Ne, Qe S (7, -y;) 5 Si, Y y, 
¡al i=l El dal 


b, b, e 
=£ DE XX ! => y, (4) 
¡=l i=l 
el estimador buscado será: 
S y, 
IE (5) 


Nx q 4) 


¡=1 , 


en la obtención de f, se está aceptando que 
el sesgo se puede corregir exclusivamente 
ajustando b,, pero ello no conduce a unos 
coeficientes b, poblacionales insesgados; en 
realidad, sólo b, es insesgado para la muestra 
utilizada. Este procedimiento de corrección 
por sesgo no altera el valor del coeficiente de 
determinación. Los otros dos procedimientos 
de ajuste de b, expuestos por McCuen et al. 
(1990) se basan en los residuos en el dominio 
logarítmico. 


Ecuaciones empíricas por probar 


Las ecuaciones de regresión potencial que han 
sido probadas y que conducen a resultados 
satisfactorios para estimar crecientes de 
diseño (Q,,) han utilizado el área de cuenca (4) 
como variable fisiográfica fundamental con 
la que están relacionados todos los procesos 
hidrológicos y a partir de la cual toman 
dimensión. Después del tamaño de la cuenca, 
la longitud del cauce principal (Lcp) también 
se ha utilizado para dar dimensión con una 
relación inversa, ya que a mayor longitud, las 
crecientes son menores (McCuen et al., 1990; 
Tasker et al., 1996; GREHYS, 1996b; Pandey 
y Nguyen, 1999). El uso de ambas variables 
origina un problema de multicolinealidad, ya 
que es común que exista correlación entre ellas 
(Montgomery et al., 2002). 

Con respecto a las variables climáticas o 
meteorológicas empleadas, por lo general se 


han considerado las predicciones asociadas 
con la lluvia máxima diaria anual (PMD,,), 
por la mayor cobertura de las estaciones 
pluviométricas y la facilidad de obtención 
de tales datos, comparados, por ejemplo, con 
los registros pluviográficos. Con base en lo 
anterior, las cuatro ecuaciones que se propone 
investigar son: 


db, 
Qr, =b,"A (6) 
Qr, =bp- A% -Lep”” (7) 
b b, 
Qr, =b,* A” -PMDF (8) 
Qy, =D, A” :Lep”? - PMDE. (9) 


La comparación entre las dos primeras 
ecuaciones busca establecer si existe mejoría 
estadística al usar dos variables predictivas; 
en cambio, la comparación entre las siguientes 
dos ecuaciones potenciales propuestas se 
orienta a definir cuál resulta más conveniente 
en asociación con la tercera variable predictiva 
involucrada. 


Índices de desempeño de las ecuaciones 
empíricas 


Cuando se tienen varias ecuaciones de RLM, 
la selección de la más conveniente se puede 
buscar con base en varios indicadores en los 
residuos (Montgomery et al., 2002); es decir, 
en las diferencias entre la variable observada 
(y) y la estimada () con la ecuación de RLM 
que se contrasta. Cuando se trabaja con los 
logaritmos de las variables, tales indicadores, 
como el coeficiente de determinación (R?) o la 
estadística C, de Mallows no son confiables, 
pues los residuos se evalúan en el dominio 
logarítmico, como ya se ha indicado. 

Por lo anterior, todas las ecuaciones 
potenciales empíricas que se obtengan a través 
de la RLM se evaluarán por medio de los 
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siguientes tres índices de desempeño promedio, 
calculados en el dominio real (Pandey y 
Nguyen, 1999): 


1. DAM es la desviación absoluta media con 
unidades de m?/s: 


1 
n—-npa 


DAM = (10) 


n 
NO, => Qr, 


i=L. 


en la cual, Q', y O son las predicciones del 
gasto observado y estimado con la ecuación 
potencial empírica; n, el número de datos 
utilizados de cada variable en este estudio (23) 
(ver cuadro 1), y npa es el número de parámetro 
de ajuste, que variará de 2a 4. 


2. EEMesel error estándar medio, con unidades 
de m?/s: 


3. EREM es el error relativo estándar medio sin 
unidades: 


(12) 


índices de 
desempeño puntual, que son los errores 


También se evaluarán tres 


relativos máximo (ERmáx), mínimo (ERmín) 
y mediano (ERmed). La expresión del error 
relativo es: 


(95, -Q;,) 


: (13) 
Or, 


ER, = 


1 


Cuando el ER resulta negativo indica 
que la estimación (0;,) resultó menor que la 
predicción observada o creciente de diseño 
calculada (Q;,, cuadro 2). Lógicamente, cuando 
ER es positivo ocurrió lo contrario. 


Información hidrométrica procesada 


En el cuadro 1 se citan las características 
generales de las 22 estaciones hidrométricas 
procesadas, cuyos datos de gasto máximo 
anual fueron obtenidos del sistema BANDAS 
(IMTA, 2002). Las estaciones están expuestas 
en orden decreciente de magnitud de cuenca. 
Los datos relativos a la estación San Francisco 
proceden del Boletín Hidrológico núm. 36 
(SRH, 1970, 1975). Se observa que los registros 
anuales de tales estaciones fluctúan de 19 
a 56 datos, con una media de 36 años. Las 
propiedades fisiográficas de área de cuenca (4) 
y longitud del cauce principal (Lcp) expuestas 
en el cuadro 1 proceden del Apéndice A de 
Escalante-Sandoval y Reyes-Chávez (2002). 
En la parte final del cuadro 1 se exponen otras 
cinco estaciones hidrométricas que serán 
empleadas en el contraste; sus registros no son 
confiables y se han integrado con base en el 
sistema BANDAS y los valores expuestos por 
Escalante-Sandoval y Reyes-Chávez (2002). En 
la figura 1 se muestran las cuencas de las 29 
estaciones hidrométricas que serán utilizadas 
dentro de la Región Hidrológica Núm. 10 
(Sinaloa). 


Predicciones del gasto máximo anual (Q,,) 


Los seis periodos de retorno que tendrán las 
crecientes de diseño que se estimarán con las 
ecuaciones empíricas por desarrollar son de 
5, 10, 25, 50, 100 y 500 años. Estos intervalos 
de recurrencia cubren la mayoría de los que 
son empleados en los dimensionamientos y 
revisiones hidrológicas de las diversas obras 
hidráulicas. Se probó la calidad estadística de 
los 23 registros procesados, verificando que 
estuvieran exentos de componentes determi- 
nísticas, para buscar obtener predicciones 
confiables. En el cuadro 2 se exponen las 
predicciones del gasto máximo anual (Q,,), 
obtenidas con base en el ajuste de la distribución 
General de Valores Extremos (GVE), a través 
de los cuatro métodos siguientes (Campos- 
Aranda, 2001): momentos, sextiles, máxima 
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Cuadro 1. Características generales de las 29 estaciones hidrométricas procesadas en la Región Hidrológica Núm. 10 (Sinaloa). 


O Roo Clave Longitud | Latitud A* Lep** aos Eo 

BANDAS wG N (km?) (km) (disponibles) 

1 Huites R. Fuerte 10037 108* 22 26" 54 26 057 267 1942-1992 (51) 
2  |San Francisco |R. Fuerte = 108* 07 2654 75311 246 1941-1973 (33) 
3 San Ignacio  |R. San Miguel 10094 107* 51 26" 53' 10 920 209 1967-1985 (19) 
4  |Santa Cruz R. San Lorenzo 10040 106* 57 247 28' 8919 180 1944-2002 (52) 
5 Guatenipa II |R. Humaya 10112 107* 13 25 21" 8 252 178 1969-2001 (32) 
6  |Jaina R. Sinaloa 10036 108* 01 25 54 8179 177 1942-1998 (56) 
7  |Palo Dulce R. Chinipas 10077 108* 26' 27" 02" 6 439 184 1958-1986 (21) 
8  |Ixpalino R. Piaxtla 10065 106" 36' DON 6 166 174 1953-1999 (45) 
9 La Huerta R. Humaya 10113 106" 42 25" 22" 6149 117 1970-1999 (28) 
10  |Toahayana R. Petatlán 10110 107" 42 Se MU 5281 131 1958-1987 (27) 
11 |Chinipas R. Oteros 10064 108* 32 2725" 5 098 131 1965-2002 (24) 
12 | Tamazula R. Tamzula 10087 106* 58' 249 56' 2241 58 1963-1999 (32) 
13 |Naranjo A. Ocoroni 10029 108* 28' 25 48" 2 064 107 1939-1984 (45) 
14 | Acatitán R. Elota 10070 106" 39 247 02” 1 884 115 1955-2002 (43) 
15 |Guamúchil R. Mocorito 10031 108* 05 25" 28" 1645 65 1940-1971 (32) 
16 |Choix R. Choix 10066 108* 19 26" 44" 1 403 82 1956-2002 (38) 
17 |Badiraguato |R. Badiraguato 10079 107* 32 25 20" 1018 59 1974-1999 (26) 
18 |El Quelite R. Quelite 10083 106* 30' 23" 30" 835 42 1961-2001 (33) 
19 | Zopilote A. Cabrera 10034 108* 21” 257 44” 666 68 1939-2001 (56) 
20 |El Bledal A. El Bledal 10027 107" 08' 247 48” 371 31 1938-1994 (56) 
21 |Pericos A. Pericos 10086 107" 41 25 07' 270 31 1961-1992 (30) 
22 |La Tina A. Sibajahui 10078 108* 37" ZO 254 30 1960-1983 (24) 
23 |Bamícori A. Barotén 10057 108* 29 26" 22" 223 20 1951-1983 (33) 
I Piaxtla R. Piaxtla = 106" 25' POS DIS07 142 1958-1973 (16) 
II |Urique II R. Urique 10100 107" 50' 27" 20" 4 000 127 1968-2002 (29) 
II | Tecusiapa R. Petatlán LS 10/22S ZOO US 90 1959-1981 (23) 
IV |Cazanate A. Álamos 10098 108" 46' 26" 36" 1813 93 1968-1985 (19) 
V  |Los Molinos |A. Los Molinos = 072 25% 42" 501 89 1958-1970 (13) 
VI |Chico Ruiz A. Mezcales 10090 107" 45 25 09" 391 25 1977-2002 (19) 


* Área de cuenca. 


** Longitud del cauce principal. 


*** Gasto máximo anual. 


verosimilitud y momentos L, adoptando los 
resultados del que condujo al menor error 
estándar de ajuste (Kite, 1977). 


Predicciones de la lluvia máxima diaria 
(PMD,,) 


Campos-Aranda (2008a) procesó los registros 
de precipitación máxima diaria anual de 
67 estaciones pluviométricas de la Región 
Hidrológica Núm. 10 (Sinaloa) y expuso sus 
predicciones de periodos de retorno de 2, 5, 10 y 


100 años. Como no existe correspondencia entre 
estos periodos de retorno y los seleccionados 
para las predicciones de gasto máximo anual, 
se decidió emplear la lluvia máxima diaria 
de 10 años en los tres primeros intervalos de 
recurrencia del gasto y la lluvia con frecuencia 
de 100 años en los siguientes tres (ecuaciones 
(8) y (9)). En cada una de las 29 cuencas de las 
estaciones hidrométricas que fueron empleadas 
se localizaron las estaciones pluviométricas 
procesadas y se buscó obtener su predicción 
promedio. Lógicamente, en las cuencas pe- 
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Figura 1. Localización en la Región Hidrológica Núm. 10 de México de las cuencas de las 29 estaciones 


hidrométricas procesadas. 


queñas sólo existe la estación pluviométrica 
ubicada en donde está la hidrométrica, es 
decir, a la salida de la cuenca, tal fue el caso 
de Bamícori, La Tina, Pericos y Choix. En 
otros casos se localizaron dos estaciones 
pluviométricas dentro de la cuenca, como 
en Badiraguato, Naranjo y Guamúchil. Por 
último, en las cuencas medianas y grandes se 
ubicaron tres o más estaciones pluviométricas. 
Los valores calculados para cada cuenca se 
tienen en las dos últimas columnas del cuadro 
Ze 


Ecuaciones empíricas obtenidas por MCR 


En los cuadros 3 y 4 se exponen los parámetros 
de ajuste obtenidos en el dominio logarítmico 
mediante MCR, para las ecuaciones (6) a (9) 
y cada uno de los seis periodos de retorno 
analizados; por ello, el primer indicador de 


desempeño es el coeficiente de determinación, 
indicado como log R”. También en tales cuadros 
se citan los otros índices de desempeño 
promedio y puntual (ecuaciones (10) a (13)). Las 
predicciones observadas y estimadas, así como 
los errores relativos obtenidos en cada una de 
las 23 estaciones hidrométricas procesadas en 
los periodos de retorno extremos analizados, es 
decir, 5 y 500 años, se muestran en los cuadros 
5 y 6. 


Ecuaciones empíricas obtenidas por MCR con 
corrección por sesgo 


Los resultados de la aplicación de la ecuación 
(5) para estimar el coeficiente b, corregido, es 
decir f, en las cuatro ecuaciones empíricas 
propuestas y para los seis periodos de retorno 
analizados, se muestran en los cuadros 7 y 8, 
en los cuales se ha suprimido log R?, por ser 
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Cuadro 2. Predicciones del gasto anual (Q,,) en m?/s en las 29 estaciones hidrométricas indicadas de la Región Hidrológica 
Núm. 10 y predicciones de la precipitación máxima diaria (PMD) en milímetros en sus cuencas respectivas. 


E EEA* Periodos de retorno (Tr) en años del Q,, Tr de la PMD 
Núm. Nombre 

(m/s) 5 10 25 50 100 500 10 100 

1 Huites 1021.6 4134 6166 10 030 14 266 20 144 44 215 95 140 
2 San Francisco 376.6 2D 3258 4925 6593 8732 16 355 80 125 
3 San Ignacio 256.4 2 065 2 607 3435 4173 5 030 7605 102 150 
4 Santa Cruz 406.0 1 434 2104 3 258 4410 5883 11 118 189 2 
5 Guatenipa II 246.9 2 466 3247 4 285 5 093 5930 8 001 121 179 
6 Jaina Sil 1 347 2029 3309 4 692 6589 14 199 117 183 
7 Palo Dulce 895.5 1274 1913 3 252 4 858 7267 18 597 110 157 
8 Ixpalino 346.4 1597 2201 S899 4 485 5 982 11 455 143 298 
9 La Huerta 96.3 1 431 1718 2 027 2 223 2 393 2 709 110 186 
10 Toahayana 158.1 17 2 019 2.653 3164 3709 9125 98 148 
11 Chinipas 139.0 1231 1626 2 209 2714 3 283 4 919 110 157 
12 Tamazula 145.0 770 1 009 1390 1745 AS 3543 145 235 
13 Naranjo 156.3 909 1 353 2 105 2 842 3775 7009 107 126 
14 | Acatitán 226.9 11227 1775) 2 654 3459 4 422 7 458 150 309 
15 Guamúchil 234.9 948 1351 2 038 2717 3579 6599 108 149 
16 Choix 99.9 461 650 979 Ta 1741 3294 98 150 
17  |Badiraguato 1077.2 1374 2 283 4 244 6664 10 381 28 652 177 371 
18 El Quelite 86.7 707 987 1414 1755 2 238 3568 150 270 
19 Zopilote 48.5 558 738 977 1161 1351 1817 121 168 
20 El Bledal 64.0 401 565 834 1090 1 404 2 446 138 247 
21 Pericos 28.5 354 457 601 719 846 1184 136 271 
22 La Tina 86.2 144 pas) 393 573 824 1 860 128 231 
23 Bamícori 52.3 270 391 591 782 1020 1815 110 202 
I Piaxtla 909.8 1672 PIS) 4 442 6613 9796 24 180 143 298 
TI Urique Il 35.2 456 549 660 740 817 983 ZO 121 
TI Tecusiapa 173.6 1 542 2 083 2861 SNS 4240 6241 97 139 
IV |Cazanate 419.6 722 1123 1907 2792 4 049 9 427 105 125 
V  |Los Molinos 12.6 226 260 293 gio 328 355 105 156 
VI  |Chico Ruiz 27:17 329 413 510 576 637 761 136 271 

* Error estándar de ajuste. 
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igual al indicado en los cuadros 3 y 4. También 
se tienen en los cuadros 7 y 8 los indicadores de 
desempeño definidos en las ecuaciones (10) a 
(13). Nuevamente, las predicciones observadas 
y estimadas, así como los errores relativos 
obtenidos en cada una de las 23 estaciones 
hidrométricas procesadas en los periodos de 
retorno extremos analizados se muestran en los 
cuadros 9 y 10. 


Ecuaciones empíricas obtenidas con 
optimización numérica 


Desde el inicio de los años noventa (McCuen 
et al., 1990; GREHYS, 1996a) se ha sugerido 


emplear la optimización numérica para 


ajustar modelos de regresión potencial 
(ecuación (1)) en el dominio real y así evitar 
los inconvenientes de la transformación 
logarítmica. También desde tal fecha se ha 
sugerido utilizar como función objetivo (FO) 
minimizar la suma de los errores relativos 
al cuadrado por medio del algoritmo de 
múltiples variables no restringidas de 
Rosenbrock (Rosenbrock, 1960; Kuester y 
Mize, 1973; Campos-Aranda, 2003). 

Con el propósito de tomar en cuenta el 
número de parámetros de ajuste (npa), el cual 
cambia en las ecuaciones (6) a (9) de dos a cua- 


tro, se consideró conveniente emplear como FO 
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Cuadro 3. Coeficientes de ajuste e índices de desempeño de las ecuaciones empíricas indicadas ajustadas por MCR 


(dominio logarítmico) para las crecientes de diseño (Q,,) en la Región Hidrológica Núm. 10. 


Coeficientes de ajuste e Periodos de retorno (Tr) en años del gasto (Q,,) 
índices de desempeño 5 10 25 50 100 500 

Q, =b, 4% 

b, 17.8007 26.7186 41.1347 54.4500 70.4608 1P22585 
b, 0.5121 0.5036 0.4981 0.4965 0.4962 0.4973 
log R? 0.8437 0.8300 0.7927 0.7487 0.6936 0.5460 
DAM 273.8 415.8 704.5 1076.5 1 681.4 4 674.3 
EEM 365.2 587.2 1 118.8 1 817.2 AUDI 8 346.5 
EREM 0.3485 0.3237 0.3270 0.3644 0.4339 0.7150 
ERmáx 1.1066 0.8969 0.6507 0.8629 1.2336 2.4565 
ERmín 0.5505 0.6172 0.6948 0.7455 0.7891 0.8664 
ERmed 0.0023 0.0082 0.0732 0.0203 0.0295 0.0087 
Q,, =b, * A*: Lep'2 

b, 17.5616 25.3697 37.3631 47.9470 60.2210 HITS 
b, 0.4974 0.4474 0.3938 0.3586 0.3259 0.2552 
d, 0.0283 0.1084 0.2012 0.2660 0.3285 0.4670 
log R? 0.8438 0.8308 0.7952 0.7528 0.6995 0.5554 
DAM 287.0 432.5 EDS 1 103.1 1 723.6 4 833.8 
EEM 375.6 609.1 1 166.0 1 894.8 3 037.6 8 664.5 
EREM 0.3578 0.3327 0.3318 0.3637 0.4266 0.6877 
ERmáx 1.1094 0.9078 0.6683 0.7483 1.0653 2.0921 
ERmín 0.5505 0.6168 0.6942 0.7449 0.7884 0.8658 
ERmed 0.0019 0.0005 0.0772 0.0244 0.0165 0.0532 


el error relativo cuadrado medio (ERCM), cuya 
expresión es: 


«(6 -05,)] 


2 


i 
¡=1 Or, 


ERCM=-" (14) 
(n-npa) 


Las condiciones iniciales que se utilizaron 
para poner a funcionar el algoritmo de 
Rosenbrock fueron las dos siguientes de 
acuerdo con cada ecuación (6) a (9), cuya FO 
se estaba minimizando: (1) se empleó el b, y los 
respectivos b, y (2) se cambió sólo al b, por f,. En 
los cuadros 11 y 12 se tienen los nuevos valores 
de los coeficientes de ajuste, así como otras 
magnitudes asociadas con el funcionamiento 
del algoritmo de Rosenbrock, como son 


FO inicial y final, número de etapas y de 
evaluaciones de la FO. Cuando se trabajó con 
las ecuaciones (6) y (7), la mayoría de las veces 
los coeficientes de ajuste de MCR condujeron a 
la FO más baja (cuadro 11). En cambio, cuando 
procesan las ecuaciones (8) y (9), en la mitad de 
los casos la segunda condición de inicio es la 
que lleva a la menor FO (cuadro 12). 


Análisis y selección de ecuaciones 
empíricas 


La comparación del desempeño entre las ecua- 
ciones (6) y (7) se tiene en el cuadro 3, e indica 
que en el dominio logarítmico, el coeficiente 
de determinación (R?) es escasamente superior 
en la ecuación (7). Los índices de desempeño 
DAM y EEM son también ligeramente mayores 
en la ecuación (7), quizás debido al mayor 
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Cuadro 4. Coeficientes de ajuste e índices de desempeño de las ecuaciones empíricas indicadas ajustadas por MCR 


(dominio logarítmico) para las crecientes de diseño (Q,,) en la Región Hidrológica Núm. 10. 


Coeficientes de ajuste e Periodos de retorno (Tr) en años del gasto (Q,,) 

índices de desempeño 5 10 25 50 100 500 
Q,, =D, * A": PMD'> PMD,, PMD 

b, 0.2490 0.2335 0.1966 1.2051 1.1086 0.8310 
b, 0.5618 0.5588 0.5604 0.5616 0.5671 0.5825 
b, 0.8116 0.9011 1.0158 0.6271 0.6832 0.8213 
log R? 0.8722 0.8657 0.8370 0.7923 0.7417 0.6004 
DAM 313.6 474.1 787.5 1 173.0 1 834.4 5 082.7 
EEM 402.8 644.4 1 209.5 1 878.5 2 999.2 8517.4 
EREM 0.3372 0.3161 0.3279 0.3749 0.4517 0.7505 
ERmáx 0.9911 0.7827 0.6737 0.9279 1.3174 2.6118 
ERmín 0.4079 0.4798 0.5686 0.6388 0.6912 0.7888 
ERmed 0.0781 0.0445 0.0207 0.0061 0.0763 0.0526 
Q,, =D, A"- Lep”2- PMD'> PMD,, PMD o 

b, 0.2482 0.2271 0.1859 0.9893 0.8717 0.5936 
b, 0.5563 0.5126 0.4671 0.4065 0.3781 0.3181 
b, 0.0107 0.0889 0.1793 0.3010 0.3666 0.5130 
b, 0.8113 0.8983 1.0102 0.6359 0.6940 0.8364 
log R? 0.8722 0.8662 0.8390 0.7977 0.7491 0.6117 
DAM 9299) 499.8 830.5 1 230.9 1 936.5 5 447.8 
EEM 413.7 668.4 1 261.3 1976.6 3151.22 8 895.2 
EREM 0.3464 0.3224 0.3257 0.3647 0.4313 0.6999 
ERmáx 0.9936 0.7919 0.5550 0.7957 1.1243 2.1997 
ERmín 0.4074 0.4797 0.5686 0.6360 0.6883 0.7859 
ERmed 0.0775 0.0536 0.0105 0.0611 0.0880 0.1448 


número de parámetros de ajuste. El índice 
adimensional EREM (ecuación (12)) indica 
que la ecuación (6) tiene mejor desempeño en 
los primeros tres periodos de retorno y en los 
siguientes tres la (7). Lo mismo se puede decir 
para la comparación entre las ecuaciones (8) 
y (9) del cuadro 4, pero la ecuación (8) sólo es 
recomendable en los dos primeros periodos de 
retorno, es decir, en 5 y 10 años. En los cuadros 
3, 4,7, 8, 11 y 12, los coeficientes de ajuste de las 
ecuaciones empíricas seleccionadas se muestran 
sombreados. 

Respecto a la comparación mostrada 
en el cuadro 7 y de acuerdo con el índice 
de desempeño EREM, se concluye que la 
ecuación (6) con f, es mejor en los cuatro 
primeros periodos de retorno y la ecuación 


(7) en los dos últimos, es decir, en 100 y 500 
años. Para la comparación del cuadro 8 y 
nuevamente con respecto al índice EREM, la 
ecuación (8) con f, resulta recomendable sólo 
en los dos primeros periodos de retorno y la 
ecuación (9) en el resto. 

El algoritmo de Rosenbrock en todos los 
casos operó y minimizó la FO inicial (ecuación 
(14)) definida con base en los resultados del 
método de MCR o de MCR con corrección por 
sesgo, sus resultados (cuadros 11 y 12), es decir, 
las ecuaciones (6) a (9) ajustadas a cada uno de 
los seis periodos de retorno serán sometidas 
al mismo proceso de selección, en este caso 
adoptando la que condujo al mejor ERCM o 
menor FO. Para el cuadro 11, la ecuación (6) 
es la mejor en los tres primeros periodos de 


Campos-Aranda, Obtención de ecuaciones empíricas para estimación de crecientes de diseño en la Región Hidrológica... 


Cuadro 5. Errores relativos de periodo de retorno de cinco años estimados en las 23 estaciones hidrométricas de la Región 


Hidrológica Núm. 10. La ecuación empírica indicada se ajustó por MCR (dominio logarítmico). 


pim | name aa 
1 Huites 4134 3 249.613 0.2139 0.5505 
2 San Francisco 229). 2 653.714 0.1574 0.3101 
3 San Ignacio 2 065 2 081.665 0.0081 0.2687 
4 Santa Cruz 1 434 1 876.696 0.3087 0.2139 
5 Guatenipa II 2 466 1 803.461 0.2687 0.2108 
6 Jaina 1 347 1 795.274 0.3328 0.1670 
7 Palo Dulce 1274 1 588.301 0.2467 0.1158 
8 Ixpalino 11557 1 553.452 0.0023 0.1094 
9 La Huerta 1 431 1 551.257 0.0840 0.0815 
10 Toahayana 17 1 434.961 0.0784 0.0784 
11 Chinipas 1231 1 409.277 0.1448 0.0243 
112 Tamazula 770 925.120 0.2015 0.0023 
13 Naranjo 909 886.951 0.0243 0.0081 
14 Acatitán 11 2z7 846.459 0.3101 0.0511 
15 Guamúchil 948 789.652 0.1670 0.0840 
16 Choix 461 727.856 0.5789 0.1448 
17 Badiraguato 1374 617.596 0.5505 0.1574 
18 El Quelite 707 557.998 0.2108 0.2015 
19 Zopilote 558 496.979 0.1094 0.2467 
20 El Bledal 401 368.310 0.0815 0.3087 
21 Pericos 354 312.996 0.1158 0.3328 
22 La Tina 144 303.356 1.1066 0.5789 
23 Bamícori 270 283.794 0.0511 1.1066 


retorno y la ecuación (7), en los tres siguientes. 
Respecto al cuadro 12, la ecuación (8) es la 
más conveniente en los dos primeros periodos 
de retorno y la ecuación (9) en los cuatro 
siguientes, es decir, en 25, 50, 100 y 500 años de 
intervalo de recurrencia. 

En el cuadro 13 se exponen las 36 
ecuaciones empíricas obtenidas bajos los tres 
planteamientos seguidos de: (1) mínimos 
cuadrados de los residuos (MCR) en el dominio 
logarítmico, (2) MCR con corrección por sesgo 
y (3) optimización numérica de una función 
objetivo (FO = ERCM). 


Contrastes numéricos de las 36 
ecuaciones empíricas seleccionadas 


Se realizaron dos contrastes numéricos. El 
primero en las seis estaciones hidrométricas 
indicadas al final de los cuadros 1 y 2, las 


cuales no fueron utilizadas en la deducción de 
las ecuaciones empíricas, porque sus registros 
no son confiables. Para el segundo contraste 
se emplearon las estaciones hidrométricas con 
registro más amplio de gasto máximo anual, 
mismo que varía de 45 a 56; se eliminó Huites 
por tener un área de cuenca que está fuera de 
las expectativas de aplicación de las ecuaciones 
empíricas dentro de la Región Hidrológica 
Núm. 10. Debido a que se han seleccionado 
seis ecuaciones empíricas para cada uno de 
los seis periodos de retorno procesados, se 
optó por obtener su valor mediano para definir 
la predicción estimada y contrastarla con la 
observada. 

En el cuadro 14 se tienen los resultados del 
primer contraste. Respecto a las predicciones 
observadas en la estación Piaxtla, éstas se 
consideran sumamente elevadas, pues son 
superiores a las observadas en la estación 
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Cuadro 6. Errores relativos de periodo de retorno de 500 años estimados en las 23 estaciones hidrométricas de la Región 


Hidrológica Núm. 10. La ecuación empírica indicada se ajustó por MCR (dominio logarítmico). 
Núm. | Nombre OS Oso = 0.5936 40981 Lep)9510 PUD 2350 si e ab EE 20 mm 
Huites 44 215 16 520.40 0.6264 0.7859 
San Francisco 16 355 12 701.47 0.2234 0.6264 
San Ignacio 7.305 11 704.96 0.6023 0.4822 
Santa Cruz 11 118 13 630.15 0.2260 0.4694 
Guatenipa II 8001 11 431.49 0.4288 0.3973 
Jaina 14 199 11 578.39 0.1846 0.2802 
Palo Dulce 18 597 9 628.99 0.4822 0.2234 
Ixpalino 11 455 15 773.59 0.3770 0.1846 
La Huerta 2709 8 667.87 2.1997 0.0696 
Toahayana 5123 7 228.64 0.4110 0.0084 
Chinipas 4919 7 509.80 0.5267 0.0394 
Tamazula 3 543 5334.19 0.5056 0.1448 
Naranjo 7 009 4 224.03 0.3973 0.2090 
Acatitán 7 458 9 016.47 0.2090 0.2260 
Guamúchil 6 599 3 501.33 0.4694 0.3770 
Choix 3294 3770.87 0.1448 0.4110 
Badiraguato 28 652 6 133.80 0.7859 0.4288 
El Quelite 3568 3 708.53 0.0394 0.5056 
Zopilote 1817 2971.44 0.6354 0.5267 
El Bledal 2 446 DAI 0.0696 0.6023 
Pericos 1 184 2 222.87 0.8774 0.6354 
22 La Tina 1860 1 875.68 0.0084 0.8774 
23 Bamícori 1815 1306.51 0.2802 2.1997 


Ixpalino, que tiene un área de cuenca mayor 
un 16%. Por lo anterior, las predicciones de 
las ecuaciones empíricas están más acordes 
con los valores probables. Lo contrario ocurre 
en la estación Urique IL, cuyas predicciones 
observadas son en extremo bajas, como se 
comprueba al compararlas con las de Chinipas, 
cuenca vecina y adyacente, con un 27% más 
de área. Nuevamente, las predicciones de 
las ecuaciones empíricas se consideran más 
aproximadas a las crecientes probables. 

En las estaciones Tecusiapa y Cazanate 
todas las predicciones de las ecuaciones 
empíricas están por debajo de las observadas, 
pero dentro del orden de magnitud. Lo 
contrario, pero mucho más acusado, ocurre 
en las dos últimas estaciones contrastadas: 
Los Molinos y Chico Ruiz. Para la primera ya 
se ha demostrado que su registro de gastos 
máximos anuales tiene valores que no están 


acordes con su tamaño de cuenca (Campos- 
Aranda, 2008b) y respecto a la segunda, la 
cuenca contigua de Pericos, que tiene un 
31% menos de tamaño, reporta predicciones 
mayores, intuyendo que quizás en el registro 
corto de Chico Ruiz todavía no se han 
presentado los valores extremos grandes, pero 
sí los menores (1991 con 10 m*/s y 1999 con 17 
m/s). Respecto al cuadro 14, las predicciones 
estimadas con las ecuaciones empíricas son 
consideradas más apegadas a la realidad que 
las derivadas de sus registros, pues éstos no 
son confiables. Las aplicaciones numéricas del 
cuadro 14 ilustran los beneficios del método 
regional desarrollado. 

En relación con los resultados del cuadro 
15, se concluye de manera general una 
excelente aproximación entre las predicciones 
observadas y las estimadas, en los periodos de 
retorno menores de 100 años. Exclusivamente, 
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Cuadro 7. Coeficientes de ajuste e índices de desempeño de las ecuaciones empíricas indicadas ajustadas por MCR con 


corrección por sesgo para las crecientes de diseño (Q,,) en la Región Hidrológica Núm. 10. 


Coeficientes de ajuste e Periodos de retorno (Tr) en años del gasto (Q,,,) 

índices de desempeño 5 10 25 50 100 500 
Op = fo: An 
de 18.3055 27.7196 43.5402 59.0646 78.9490 153.4207 
b, 0.5121 0.5036 0.4981 0.4965 0.4962 0.4973 
DAM 281.0 432.0 745.3 ZO MUZEMS 5 029.8 
EEM 362.6 579.8 1095.9 1766.9 2 813.8 7910.6 
EREM 0.3638 0.3439 0.3619 0.4215 0.5277 1.0002 
ERmáx 1.1664 0.9680 0.7472 1.0208 1.5027 3.3376 
ERmín 0.5378 0.6028 0.6769 0.7240 0.7637 0.8323 
ERmed 0.0260 0.0290 0.1360 0.1068 0.1535 0.2440 
Q, =f: A% Lep” 
f, 18.0685 26.3414 39.5923 52.0804 67.5652 122.8478 
b, 0.4974 0.4474 0.3938 0.3586 0.3259 0.2552 
b, 0.0283 0.1084 0.2012 0.2660 0.3285 0.4670 
DAM 258) 472.3 788.8 ISS 1 860.7 5 247.5 
EEM 372.8 601.7 1 141.4 1 841.7 2 931.8 8 224.2 
EREM 0.3737 0.3692 0.3682 0.4227 0.5226 0.9733 
ERmáx 1.1703 1.0259 0.7678 0.8991 1.3171 2.8842 
ERmín 0.5375 0.5931 0.6760 0.7229 0.7626 0.8314 
ERmed 0.0308 0.0625 0.1415 0.1127 0.1405 0.1893 


en la estación Naranjo se tienen errores 
relativos por defecto mayores del 20% en 
los periodos de retorno mayores de 50 años. 
En el resto de las estaciones hidrométricas 
contrastadas, los errores relativos por exceso 
no son importantes, pues llegan a un máximo 
del 38% en Santa Cruz en el periodo de retorno 
de cinco años. En las estaciones Ixpalino y El 
Bledal se tienen los errores relativos más bajos 
en todas sus predicciones, en la primera por 
exceso y en la segunda por defecto, pero ambos 
menores del 8.5%. 


Conclusiones 


Se estableció encontrar sólo cuatro tipos de 
ecuaciones potenciales para las crecientes de 
diseño de periodos de retorno de 5, 10, 25, 50, 
100 y 500 años, empleando como variables 
predictivas dos propiedades fisiográficas, el 
área de cuenca y la longitud del cauce principal, 
y una característica climática, la precipitación 


máxima diaria de periodo de retorno de 10 y 
100 años. Sin embargo, sólo se seleccionan dos 
ecuaciones para cada intervalo de recurrencia: 
las 
emplear tres procedimientos de obtención de 


de mejor desempeño estadístico. Al 


los parámetros de ajuste, se obtiene un total 
de 36 ecuaciones empíricas, expuestas en el 
cuadro 13. 

Los dos primeros métodos de ajuste, el de mí- 
nimos cuadrados de los residuos en el dominio 
logarítmico y su versión corregida por sesgo 
enel dominioreal conducen a cuatro ecuaciones 
empíricas en cada periodo de retorno, que son 
únicas; por el contrario, el tercer procedimiento 
de ajuste, basado en optimización numérica 
(algoritmo de Rosenbrock) permite obtener 
muchas ecuaciones empíricas con sólo cambiar 
la función objetivo (FO) a minimizar. En este 
trabajo, únicamente se aplicó como FO el error 
relativo cuadrado medio. 

Los dos contrastes realizados, el primero 
en seis estaciones hidrométricas no usadas en 
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Cuadro 8. Coeficientes de ajuste e índices de desempeño de las ecuaciones empíricas indicadas ajustadas por MCR con 


corrección por sesgo para las crecientes de diseño (Q,,) en la Región Hidrológica Núm. 10. 


Coeficientes de ajuste e Periodos de retorno (Tr) en años del gasto (Q,,) 
índices de desempeño 5 10 25 50 100 500 
Q,=f,* A*- PMD' PMD,, PMD., 
de 0.2567 0.2424 0.2078 112975 Peso 1.0285 
b, 0.5618 0.5588 0.5604 0.5616 0.5671 0.5825 
b, 0.8116 0.9011 1.0158 0.6271 0.6832 0.8213 
DAM 308.4 476.6 802.6 1227.6 1 966.0 5 554.8 
EEM 398.4 635.1 1 185.2 1 833.0 2 905.4 8 105.0 
EREM 0.3534 0.3360 0.3611 0.4257 0.5375 1.0213 
ERmáx 1.0526 0.8507 0.7690 1.0757 1.5736 3.4702 
ERmín 0.3896 0.4600 0.5441 0.6112 0.6571 0.7386 
ERmed 0.0495 0.0080 0.0351 0.0701 0.1954 0.3028 
Q,,=f,* A'1- Lep'2- PMD' PMD,, PMD,., 
A 0.2556 0.2358 0.1966 1.0651 0.9680 0.7329 
b, 0.5563 0.5126 0.4671 0.4065 0.3781 0.3181 
b, 0.0107 0.0889 0.1793 0.3010 0.3666 0.5130 
b, 0.8113 0.8983 1.0102 0.6359 0.6940 0.8364 
DAM 324.6 502.7 846.5 1 299.4 ANDAN 5 953.0 
EEM 409.4 658.9 1 236.8 1 932.9 3 062.6 8 522.5 
EREM 0.3624 0.3429 0.3598 0.4156 0.5168 0.9620 
ERmáx 1.0530 0.8606 0.6445 0.9333 1.3589 2.9205 
ERmín 0.3897 0.4598 0.5438 0.6081 0.6539 0.7357 
ERmed 0.0500 0.0174 0.0686 0.1424 0.2082 0.4134 
la deducción de las ecuaciones empíricas y, el Referencias 


segundo, en las seis estaciones de aforos con 
mayor registro, destacan la regularidad de las 
predicciones estimadas, cuyo valor mediano 
se considera una aproximación aceptable a los 
valores probables de las crecientes de diseño en 
cualquier sitio y su cuenca dentro de la Región 
Hidrológica Núm. 10 (Sinaloa) de México. Por 
lo anterior, el procedimiento seguido para 
obtener las ecuaciones potenciales que permiten 
estimar predicciones de gasto máximo anual 
o crecientes de diseño se sugiere aplicarlo 
o desarrollarlo en otras regiones o zonas 
geográficas para obtener tales herramientas de 
estimación hidrológica (ecuaciones empíricas), 
cuya aplicación es extraordinariamente simple. 
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Cuadro 9. Errores relativos de periodo de retorno de cinco años estimados en las 23 estaciones hidrométricas de la Región 
Hidrológica Núm. 10. La ecuación empírica indicada se ajustó por MCR y se corrigió por sesgo. 


E Predicciones Error relativo ER 
as noi observadas LS (ER) ordenado 

1 Huites 4 134 3 341.767 0.1916 0.5378 
2 San Francisco 21292 2.727.941 0.1902 0.2906 
3 San Ignacio 2 065 2 140.698 0.0367 0.2479 
4 Santa Cruz 1 434 1929.917 0.3458 0.1916 
5 Guatenipa II 2 466 1 854.604 0.2479 0.1884 
6 Jaina 1 347 1 846.185 0.3706 0.1434 
7 Palo Dulce 1274 1 633.342 0.2821 0.0908 
8 Ixpalino 1557 1 597.505 0.0260 0.0841 
9 La Huerta 1 431 1 595.249 0.1148 0.0555 
10 Toahayana 1557 1 475.654 0.0522 0.0522 
11 Chinipas 1231 1 449.242 0.1773 0.0034 
162 Tamazula 770 SILOS 0.2355 0.0260 
13 Naranjo 909 912.104 0.0034 0.0367 
14 Acatitán 1227 870.463 0.2906 0.0809 
15 Guamúchil 948 812.045 0.1434 0.1148 
16 Choix 461 748.496 0.6236 0.1773 
17 Badiraguato 1374 635.110 0.5378 0.1902 
18 El Quelite 707 573.822 0.1884 0.2355 
19 Zopilote 558 511.073 0.0841 0.2821 
20 El Bledal 401 378.755 0.0555 0.3458 
21 Pericos 354 321.872 0.0908 0.3706 
22 La Tina 144 SOS 1.1664 0.6236 
23 Bamícori 270 291.842 0.0809 1.1664 
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Cuadro 10. Errores relativos de periodo de retorno de 500 años estimados en las 23 estaciones hidrométricas de la Región 
Hidrológica Núm. 10. La ecuación empírica indicada se ajustó por MCR y se corrigió por sesgo. 


Núm. | Nombre decia Qs = 0.7329 A03I Lep0510 PMD 259% dde E mis a e e 
1 Huites 44 215 20 397.250 -0.5387 0.7357 
2 San Francisco 16 355 15 682.120 0.0411 0.5387 
3 San Ignacio 7 305 14 451.770 0.9783 0.3607 
4 Santa Cruz 11 118 16 828.740 0.5136 0.3449 
5 Guatenipa II 8001 14 114.110 0.7640 -0.2559 
6 Jaina 14 199 14 295.490 0.0068 0.1112 
7 Palo Dulce 18 597 11 888.630 0.3607 0.0411 
8 Ixpalino 11 455 19 475.170 0.7001 0.0068 
9 La Huerta 2709 10 701.950 2.9505 0.1488 
10  |Toahayana Bi23 8 924.988 0.7421 0.2451 
11  |Chinipas 4919 9 272.124 0.8850 0.2833 
12 | Tamazula 3 543 6 585.969 0.8589 0.4134 
13 [Naranjo 7009 5 215.279 0.2559 0.4927 
14 | Acatitán 7 458 11 132.360 0.4927 0.5136 
15  |Guamúchil 6599 4 322.990 0.3449 0.7001 
16  |Choix 3294 4 655.785 0.4134 0.7421 
17 — |[Badiraguato 28 652 7 573.219 0.7357 0.7640 
18 [El Quelite 3568 4 578.810 0.2833 0.8589 
19 | Zopilote 1817 3 668.745 1.0191 0.8850 
20 |El Bledal 2 446 2 809.847 0.1488 0.9783 
21  |Pericos 1 184 2 744.513 1.3180 1.0191 
22  |LaTina 1860 2 315.850 0.2451 1.3180 
23  |Bamícori 1815 1 613.112 0.1112 2.9505 
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Cuadro 11. Coeficientes de ajuste y funciones objetivo (FO) de las ecuaciones empíricas indicadas ajustadas con optimización 


numérica (algoritmo de Rosenbrock) para las crecientes de diseño (Q,,) en la Región Hidrológica Núm. 10. 


Coeficientes de ajuste e 
índices de desempeño 


Periodos de retorno (Tr) en años del gasto (Q,,) 


5 10 25 50 100 500 
Q,=b,: An 
b, 9.6794 16.6157 31.8794 53.0888 66.9377 110.4147 
b, 0.5712 0.5488 0.5137 0.4779 0.4714 0.4542 
FO inicial 0.1215 0.1048 0.1070 0.1328 0.1883 0.5112 
Número de etapas 8 13 7 3 2 4 
Evaluaciones de la FO 48 70 41 22 18 30 
FO final 0.0913 0.0818 0.0848 0.1029 0.1327 0.2536 
Q;,= by: A'I: Lep” 
b, 15.1401 23.7836 35.4876 34.6689 57.1394 88.0143 
b, 0.7466 0.4465 0.3899 0.2958 0.3091 0.2327 
b, 0.3970 0.0978 0.1914 0.4117 0.3245 0.4376 
FO inicial 0.1397 0.1107 0.1101 0.1322 0.1820 0.4730 
Número de etapas 11 4 3 13 3 3 
Evaluaciones de la FO 96 43 45 117 37 33 
FO final 0.0983 0.0924 0.0901 0.1015 0.1310 0.2481 


Cuadro 12. Coeficientes de ajuste y funciones objetivo (FO) de las ecuaciones empíricas indicadas ajustadas con optimización 


numérica (algoritmo de Rosenbrock) para las crecientes de diseño (Q,,) en la Región Hidrológica Núm. 10. 


Coeficientes de ajuste e 


Periodos de retorno (Tr) en años del gasto (Q,,) 


índices de desempeño 5 10 25 50 100 500 
Q,,=b,* A*- PMD? PMD,, PMD o 
b, 0.1821 0.2303 0.1971 1.1448 1.1822 0.9922 
b, 0.6011 0.5728 0.5747 0.5335 0.5436 0.4948 
b, 0.7873 0.8560 0.9641 0.6428 0.6649 0.8342 
FO inicial 0.1137 0.1129 0.1076 0.1405 0.2889 1.0430 
Número de etapas 7 Pl 12 p) 4 7 
Evaluaciones de la FO 69 25 118 21 58 73 
FO final 0.0914 0.0826 0.0879 0.1101 0.1442 0.2751 
Q, =0b,* A: Lep??- PMD' PMD,, PMD o 
b, 0.1953 0.5055 0.1966 0.9856 0.9080 0.4005 
b, 0.5755 0.5367 0.4034 0.3911 0.3636 0.1872 
b, 0.0502 0.0889 0.3017 0.3511 0.3844 0.7050 
d, 0.7667 0.6674 0.9616 0.5876 0.6559 0.8674 
FO inicial 0.1200 0.1176 0.1294 0.1730 0.2671 0.4899 
Número de etapas 8 21 10 8 5 16 
Evaluaciones de la FO 118 25 150 7 79 180 
FO final 0.0962 0.0841 0.0853 0.1034 0.1350 0.2540 
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Cuadro 13. 36 ecuaciones empíricas obtenidas con los planteamientos indicados, para los seis periodos de retorno indicados. 


Ajuste por MCR en el dominio 
logarítmico 


Ajuste por MCR con corrección por 
sesgo 


Ajuste por optimización numérica de 
una FO 


Q, =17.8007 A052 


Q,= 18.3055 A052! 


O, = 9.6794 4037 


Q 1, = 26.7186 4030 


Q1) = 27.7196 A03%6 


Qy, =16.6157 A05%8 


Q,; = 41.1347 40051 


O, = 43.54.02 40481 


Q,, = 31.8794 40517 


es == 47.9470 403586 IL gueka 


Qs) = 59.0646 A01965 


(De = 34.6689 402958 IL pui 


Dio an 60.2210 403259 Lepo3285 


Oi =67.5652 403259 Lep032s 


Quo =57.1394 403091 Lepo225 


Qs = 97.7937 ASEO 


(e = 122.8478 A02552 IL pesto 


(OL. ES ¿0 Lg ESE 


Q, = 0.2490 APEPMD AS 


O, = 0.2567 A0%5 PMDOs6 


Q, == 0.1821 Ao0so1 PMDI 


Qs, = 0.2335 40558 PMDO9908 


Q;o = 0.2424. 40558 PMDO901 


Q io = 0.2303 49578 PMD25%0 


OQ), = 0.1859 AQ4671 Epia PMD ¡22 


O); an 0.1966 A04671 Lepo ra PMD2 


Os == 0.1966 A04034 Lepe PMD 7.9616 


(Qu = 0.9893 A0-4065 Ines LID 


om = 1.0651 04065 IGprson LID) 


Qu = 0.9856 pos IL pese LD O 


Qi0o == 0.8717 A03781 Ecpote PMD ¿29 


Quo; = 0.9680 A03781 Ecos PMD 


Dio = 0.9080 A03636 Lprtte PMD 22 


(e =; 0.5936 A03181 Jr PAIS 


a =; 0.7329 A03181 ILGpites ID 


Quo = 04005 401872 T cp07050 PM 9.5674 


Cuadro 14. Predicciones estimadas (m?/s) con las ecuaciones empíricas seleccionadas en seis estaciones hidrométricas de 


registro no confiable de la Región Hidrológica Núm. 10. 


Piaxtla 


Urique II 


Ecuaciones - a A 
70 Periodos de retorno en años Periodos de retorno en años 
empíricas del 
5 10 25 50 100 500 5 10 25 50 100 500 
Cuadro 3 1439 | 2008 | 2948 | 3882 | 5020 8831 | 1245 | 1741 | 2561 | 3405 | 4413 | 7799 
Cuadro 4 1780 | 248 | 378 | 970 | 7485 | (8530 874 1177 | 1672 | 2614 | 3304 | 5503 
Cuadro 7 1479 | 2083 | 3121 | 4176 | 5632 | 11094 | 1280 | 1806 | 2711 | 3629 | 4952 | 9797 
Cuadro 8 1784 | 2560 | 3951 | 5791 | 7950 | 16729 901 1222 | 1768 | 2814 | 3669 | 6794 
Cuadro 11 1299 | 1840 | 2612 | 3372 | 4043 5665 | 1105 | 1575 | 2259 | 2962 | 3573 | 5051 
Cuadro 12 1571 | 219 | 3297 | 4571 | 5789 9192 798 15 | 1629) | 287 | 27% 3 687 
Q,, mediano 1525 | 2138 | 3209 | 4374 | 5711 | 10143 | 1003 | 1399 | 2014 | 2888 | 3621 5923 
Q,, observado 1672 | 2585 | 4442 | 6613 | 9796 | 24180 456 549 660 740 817 983 
. Tecusiapa Cazanate 
Ecuaciones 7 A 
¿ala Periodos de retorno en años Periodos de retorno en años 
empíricas del 
5 10 25 50 100 500 5 10 25 50 100 500 
Cuadro 3 1208 | 1691 | 2488 | 3042 | 3867 | 6542 830 1169 | 1727 | 2860 | 95078 | Sii 
Cuadro 4 1042 | 1436 | 1983 | 2513 | 3136 | 5084 737 1024 | 1535 | 1761 | 2235 | 3747 
Cuadro 7 1242 | 1754 | 2633 | 3525 | 4338 | 8218 854 1213 | 1828 | 2450 | 3454 | 6921 
Cuadro 8 1075 | 1491 | 2097 | 2706 | 3483 | 6277 759 1063 | 1623 | 189 | 2482 | 4626 
Cuadro 11 1069 | 1526 | 2192 | 2526 | 3138 | 4286 703 1020 1504 [2062 [| 2529 | 3.666 
Cuadro 12 943 1293 | 1724 | 2177 | 2602 | 3227 646 910 1398 | 1554 | 1883 | 2626 
Q,, mediano wa | 150 | 2145 | 266 | Ss | Sc 748 loe [0 | 170 | 250% | 41 
Q,, observado 1542 | 2083 | 2861 | 3516 | 4240 | 6241 | 1542 | 2083 | 2861 | 3516 | 4240 | 6241 
Los Molinos Chico Ruiz 


Ecuaciones 
empíricas del 


Periodos de retorno en años 


Periodos de retorno en años 


5 10 25 50 100 | 500 5 10 25 50 100 | 500 
Cuadro 3 430 | 612 | om 11147 | 14688 | 2514 | 378 | 540 | 804 | 960 | 1213 | 2017 
Cuadro 4 358 | 49 | 706 | 8906 | 1120 | 1815 | 384. | 59 | 769 | 1040 | 1323 | 2240 
Cuadro 7 442 | 635 | 963 | 129% | 16488 | 3158 | 389 | 560 | 851 | 1144 | 1361 | 2534 
Cuadro 8 369 | 58 | 747 | 964 | 1244 | 2240 | 39068 | 570 | 814 | 1120 | 1469 | 2765 
Cuadro 11 337 | 50 | 777 | 92 | 1237 | 1772 | 293 | 40 | 684 | 763 | 1028 | 1444 

Cuadro 12 298 | 435 | 620 | 759 937 | 1256 | 315 | 471 | 350 | 87 | 1081 | 1527 
Q,, mediano 364 | su | 762 | 953 | 1241 | 2028 | 387 | 545 | 792 | 1000 | 1268 | 2129 
Q,, observado | 226 | 260 | 293 312 | 328 | 355 | 329 | 413 | 510 | 576 | 637 | 761 
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Cuadro 15. Contraste de predicciones estimadas (m*/s) con las ecuaciones empíricas seleccionadas en seis estaciones 
hidrométricas de registro confiable de la Región Hidrológica Núm. 10. 


Santa Cruz Jaina 
al Periodos de retorno en años Periodos de retorno en años 
5 10 25 50 100 500 15 10 25 50 100 500 
Cuadro 3 1877 | 2607 | 3818 | 4980 | 6427 | 11263 | 1795 | 2496 | 3657 | 4806 | 6214 | 10930 
Cuadro 4 2184 | 3087 | 4617 | 5762 | 7528 | 13630 | 1875 | 2620 | 3884 | 5025 | 6516 | 11578 
Cuadro 7 1930 | 2705 | 4042 | 5403 | 7211 | 14148 | 1846 | 2590 | 3871 | 5176 | 6972 | 13731 
Cuadro 8 2251 | 3205 | 4883 | 6204 | 8359 | 16829 | 1933 | 2720 | 4107 | 5411 | 7236 | 14296 
Cuadro 11 1747 | 2446 | 3410 | 4335 | 5126 | 7091 | 1663 | 2332 | 3262 | 4196 | 4964 | 6899 


Cuadro 12 2028 | 2774 | 4072 | 4996 | 6145 8948 | 1740 | 2365 | 3459 | 4392 | 5356 7 627 


Q,, mediano 1979 | 2740 | 4057 | 5200 | 6819 | 12447 | 1821 | 2543 | 3764 | 4916 | 6365 | 11254 


Q,, observado 1434 | 2104 2029 | 3309 | 4692 | 6589 | 14199 


Error relativo (%)| 38.0 30.2 25.3 13.8 4.8 3.4 -20.7 
Ixpalino Naranjo 
On Periodos de retorno en años Periodos de retorno en años 
5 10 10 25) 50 100 500 
Cuadro 3 1554 | 2165 1248 | 1842 | 2566 | 3363 | 6081 
Cuadro 4 1882 | 2681 1120 | 1704 | 1946 | 2485 | 4224 
Cuadro 7 1598 | 2246 1294 | 1950 | 2613 | 3773 | 7638 
Cuadro 8 1941 | 2783 116 | 1302 | 20% [| 269 [| 25 
Cuadro 11 1415 | 1997 109 | 1608 | 2270 | 2755 | 4017 


Cuadro 12 ¡ISA (2389 996 1565 | 1725 | 2094 | 2990 


Q,, mediano 1659 | 2318 1142 | 1753 | 2183 | 2757 | 4720 


Q,, observado 57 | 22401 1555 | 210 [| 282 155 | 7005 


Error relativo (%)| 6.6 5.3 15.6 16.7 23.2 27.0 32.7 


Zopilote El Bledal 
a Periodos de retorno en años Periodos de retorno en años 

5 10 10 25 50 100 500 

Cuadro 3 497 706 526 784 997 | 1279 | 2200 
Cuadro 4 471 665 540 TED 1024 | 1316 | 2276 z 
Cuadro 7 54 732 545 829 1114 | 1436 | 2764 2 
Cuadro 8 485 690 1109 | 1386 | 1860 | 3669 389 561 837 1102 | 1461 | 2810 E 
Cuadro 11 397 589 813 | 1348 | 1676 | 2532 284 427 666 820 | 1084 | 1567 E 
Cuadro 12 396 579 973 1119 | 1408 | 2256 309 463 688 848 | 1084 | 1623 B 
Q,, mediano 478 678 1049 | 1367 | 1768 | 3320 373 533 788 1011 | 1298 | 2238 E 
Q,, observado 558 738 77 | MIG | ASS | LE 401 565 834 | 109 | 1404 | 2446 5 
Error relativo (%)| -14.3 | 8.1 74 17.7 30.9 82.7 7.0 5.7 5.5 7.2 7.5 8.5 a 
zo 
3 


Tecnología y 


